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Abstract 

We classify the unitary, renormahzable, Lorentz violating quantum field theories of interacting scalars 
and fermions, obtained improving the behavior of Feynman diagrams by means of higher space derivatives. 
Higher time derivatives are not generated by renormalization. Renormalizability is ensured by a "weighted 
power counting" criterion. The theories contain a dimensionful parameter A^, yet a set of models are 
classically invariant under a weighted scale transformation, which is anomalous at the quantum level. 
Formulas for the weighted trace anomaly are derived. The renormalization-group properties are studied. 
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1 Introduction 



The set of power-counting renormalizable theories is considerably restricted by the assumptions of 
unitarity, locahty, causahty and Lorentz invariance. If we relax one or some of these assumptions 
we can enlarge the set of renormalizable theories. However, usually the enlargement is too wide. 
For example, there exist an infinite set of renormalizable nonunitary theories. Improving the 
behavior of propagators at large momenta with the help of higher-derivative kinetic terms [T] it is 
possible to define a renormalizable higher-derivative version of every theory, including gravity p]. 
Relaxing locality can in principle make every theory renormalizable, smoothing away the small 
distance singularities that originate the UV divergences |3j. Unitarity violations due to higher 
derivatives can in some cases be traded for causality violations [UlS]. 

The purpose of this paper is to investigate the issue of renormalizability in the presence of 
Lorentz violations, while preserving both locality and unitarity. The UV behavior of propagators 
is improved with the help of higher space derivatives. It is proved that, under certain conditions, 
renormalization does not turn on terms with higher time derivatives, thus preserving unitarity. 
Renormalizability follows from a modified power-counting criterion, which weights time and space 
differently. The set of consistent theories is still very restricted, yet considerably larger than the 
set of Lorentz invariant theories. Renormalizable models exist in arbitrary spacetime dimensions. 

The quadratic terms that contain higher space derivatives, as well as certain vertices, are 
multiplied by inverse powers of a scale A/,. Despite the presence of the dimensionful parameter 
A^, certain models have a weighted scale invariance, which is anomalous at the quantum level. 
The weighted trace anomaly is worked out explicitly. 

In this paper we concentrate on scalar and fermion theories, leaving the study of gauge theo- 
ries and gravity to separate publications. Lorentz violating models with higher space derivatives 
might be useful to define the ultraviolet limit of theories that are otherwise nonrenormalizable, 
including quantum gravity, and allow to remove the divergences with a finite number of indepen- 
dent couplings. Other domains where the models of this paper might find applications are Lorentz 
violating extensions of the Standard Model [6], effective field theory [7], renormalization-group 
(RG) methods for the search of asymptotically safe fixed points [8], nonrelativistic quantum field 
theory for nuclear physics [9], condensed matter physics and the theory of critical phenomena 
|10j . Certain (/9^-models that fall in our class of renormalizable theories are useful to describe the 
critical behavior at Lifshitz points [H] and have been widely studied in that context p2], with a 
variety of applications to real physical systems. Effects of Lorentz and CPT violations on stabil- 
ity and microcausality have been studied [13], as well as the induction of Lorentz violations by 
the radiative corrections |14| . The renormalization of gauge theories containing Lorentz violating 
terms has been studied in [15]. For a recent review on astrophysical constraints on the Lorentz 
violation at high energy see ref. |16| . 



2 



The paper is organized as follows. In section 2 we study the renormalizability of scalar the- 
ories, while in section 3 we include the fermions. In section 4 we analyze the divergent parts 
of Feynman diagrams and their subtractions. We prove the locality of counterterms and study 
the renormalization algorithm to all orders. The one-loop divergences are computed explicitly. 
In section 5 we analyze the renormalization structure and the renormalization group. In section 
6 we study the energy-momentum tensor, the weighted scale invariance and the weighted trace 
anomaly. In section 7 we generalize our results to nonrelativistic theories. Section 8 contains 
the conclusions. In the appendices we collect more observations about the cancellation of sub- 
divergences and the locality of counterterms, and some expressions of Euclidean propagators in 
coordinate space. 

Preliminaries. We use the dimensional-regularization technique whenever possible. Since the 
analysis of divergences is the same in the Euclidean and Minkowskian frameworks, we write our 
formulas directly in the Euclidean framework, which is more explicit. Yet, with an abuse of 
language, we still speak of "Lorentz symmetry", since no confusion is expected to arise. 

We first consider models where the ci-dimensional spacetime manifold is split into the 
product M^(gi M-^ of two submanifolds, a d-dimensional submanifold Mj, containing time and 
possibly some space coordinates, and a d-dimensional space submanifold M-^. Lorentz and rota- 
tional symmetries in the two submanifolds are assumed. This kind of splitting could be useful 
to describe specific physical situations (for example the presence of a non isotropic medium in 
condensed matter physics), but here it is mainly used as a starting point to illustrate our argu- 
ments in concrete examples. Indeed, most Lorentz violating theories contain a huge number of 
independent vertices, so it is convenient to begin with models where unnecessary complicacies 
are reduced to a minimum. The extension of our construction to the most general case, which is 
rather simple, will be described later. In the same spirit, a number of discrete symmetries, such 
as parity, time reversal, ip —ip, etc., are often assumed. 

To apply the dimensional-regularization technique, both submanifolds have to be continued 
independently. The total continued spacetime manifold Md is therefore split into the product 
Mq (8) Mjy, where D = d — ei and D = d~ 82 are complex and D = D + D. Each momentum p is 
split into "first" components p, which live in M^, and "second" components p, which live in M^: 
p = {p,p). The spacetime index /Li is split into hatted and barred indices: fi = Notations 
such as pp, p^ and p^ refer to the same object, as well as p^, p^, pp. Frequently, Latin letters are 
used for the indices of the barred components of momenta. Finally, A = didi. 

We say that Pk,n{p,p) is a weighted polynomial in p and p, of degree k and weight 1/n, where 
A; is a multiple of 1/n, if -Pfe,n(C"P) Cp) is a polynomial of degree kn in ^. Clearly, 

Pki,n{P^P)Pk2,n{P^P) = Pki+k2,n{P^P)- 
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We say that Hk^n{p,p) is a homogeneous weighted polynomial in p and p, of degree k and weight 
1/n, if Hi^,n{\p^\^/'^p) = X''Hk^n{p,p)- It is straightforward to prove that a weighted polynomial 
Pk,n of degree k can be expressed as a linear combination of homogeneous weighted polynomials 
Hk'^n of degrees k' < k. 



2 Renormalizability by weighted power counting 

In this section we classify the renormalizable Lorentz violating scalar field theories that can be 
constructed with the help of quadratic terms containing higher space derivatives and prove that 
renormalization does not generate higher time derivatives. 

Consider a generic scalar field theory with a propagator defined by the quadratic terms 

c,,,^ = ^0^f + ^^ir^f, (2.1) 

where is an energy scale. Up to total derivatives it is not necessary to specify how the 2n 
derivatives d contract among themselves. The n of (12. Ih should be understood as the highest 
power of d that appears in the quadratic terms of the total lagrangian. Other quadratic terms of 
the form 

(B""^)^, m<n, (2.2) 



could be present, or generated by renormalization. They are weighted monomials of degrees 
< 2 and weight 1/n. For the purposes of renormalization, it is convenient to consider such 
terms as "interactions" (two-leg vertices) and treat them perturbatively. Indeed, the counterterms 
depend polynomially on the parameters Om, because when the integral associated with a graph is 
differentiated a sufficient number of times with respect to the a^'s it becomes overall convergent. 
The Om-polynomiality of counterterms generalizes the usual polynomiality in the masses. Thus 
we can assume that the propagator is defined by (|2.ip and treat every other term as a vertex. 
Then the propagator is the inverse of a weighted homogeneous polynomial of degree 2 and weight 
1/n. The coefficient of the term {d^(p)'^ must be positive, to have an action bounded from below 
in the Euclidean framework or, equivalently, an energy bounded from below in the Minkowskian 
framework. 

Label the vertices that have (/3-legs with indices a, to distinguish different derivative struc- 
tures. Each vertex of type {N, a) defines a monomial in the momenta of the fields. Denote the 
weighted degree of such a monomial by A vertex with pi derivatives d, p2 derivatives d and 



A'" 99-legs is symbolically written as 
and its weighted degree is 



Sn =Pi + —. 

n 
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Consider a Feynman graph G made of L loops, E external legs, / internal legs and vertices 
of type (A^, a). The integral associated with G has the form 

where p are the loop momenta, k are the external momenta, „ are the propagators, which 
have weighted degree —2, and V^/\^^ are the vertices, with weighted degrees Sj. The integral 

measure d^p d^p is a weighted measure of degree D= D + D/n. Performing a rescahng 
{k,k) — > {\k, X^^^k), accompanied by an analogous change of variables {p,p) (Xp, X^^"'p), 
it is straightforward to prove that Toik) is a weighted function of degree 

V 



LB-2I + ^5j = LB-2I+ ^ 5 



N '"n 



By the locality of counterterms, once the subdivergences of G have been inductively subtracted 
away, the overall divergent part of G is a weighted polynomial of degree 

u{G) = Ld-2I+ ^ 4)4") 

{N,a) 

in the external momenta, where d= d + d/n. The usual relations 

L = I-V+1, E + 2I=J2NvP, (2.3) 



(JV,a) 



allow us to write 



where 



LuiG) = d{E) + [4"^ - d{N)\ , (2.4) 

{N,a) 

d{X) = d(^l-^^+X ; (2.5) 

The theory is i) renormahzable, if it contains all vertices with 4^ < d{N), and only those: 
Lo{G) does not increase when the number of vertices increases; ii) super-renormalizable, if it 
contains all vertices with < d{N), and only those: uj{G) decreases when the number of 
vertices increases; in) strictly-renormalizable, if it contains all vertices with 4^ = d{N), and 
only those: oj{G) does not depend on v^^^; iv) nonrenormalizable, if it contains some vertices 
with 4^ > d{N): u{G) increases when the number of those vertices increases. 

The vertices with = d{N) are called "weighted marg inal", those with < d{N) are 
called "weighted relevant" and those with 4^ > d{N) are called "weighted irrelevant". 
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By locality, 6^^^ cannot be negative. Moreover, polynomiality demands that there must exist 
a bound A^max on the number of legs that the vertices can contain. It is easy to show that these 
requirements are fulfilled if and only if 

d>2 (2.6) 

and the bound is 

2d 



-^max 



where [x] denotes the integral part of x. The existence of nontrivial interactions (A'^max > 3) 
requires d < 6, while the existence of nontrivial even interactions (A'max > 4) requires d < 4. 

To complete the proof of renormalizability, observe that when < d{N) the weighted degree 
of divergence u!{G) of a graph G satisfies 

Lo{G) < d{E). (2.8) 

The inequality (12. 6p ensures also that w(G) decreases when the number of external legs increases. 
Finally, since the vertices that subtract the overall divergences of G are of type (£', q) with = 
a;(G), it is straightforward to check that the lagrangian contains all needed vertices. Indeed, (|2.8I) 

(a) 

coincides with the inequality satisfied by (5^ . 

Now we prove that the renormalizable models just constructed are perturbatively unitary, in 
particular that no higher time derivatives are present, both in the kinetic part and in the vertices, 
and no higher time derivatives are generated by renormalization. Indeed, a lagrangian term with 
higher time derivatives would have (jj^^ > 2 for > 2 or > 2 (terms with = 1 need not be 
considered, since they cannot contain derivatives). This cannot happen in a renormalizable theory, 
because (|2.6I1 and (^j^^ < d{N) imply 5^^^ < 2 in general and 6^^^ < 2 for A^ > 2. In particular, 
true vertices (A^ > 2) cannot contain any 5-derivative at all, because invariance under the reduced 
Lorentz and rotational symmetries of Mq and Mjj exclude also terms containing an odd number 
of d's or an odd number of 5's. Similar conclusions apply to the counterterms, because of (12. Sp . 
Therefore, renormalization does not turn on higher time derivatives, as promised. 

Weighted scale invariance. The strictly renormalizable models have the 5^^^ = d{N). Their 
lagrangian has the form 



{d,d) 2^ 2A2"~2 ^ {n-l){N+d-dN/2) I ^ \a ^ ' 



Here 



{N,a) NlAl 

Qnd(N)^j^ denotes a basis of lagrangian terms constructed with N fields (p and nd{N) 



9-derivatives acting on them, contracted in all independent ways, and \(^N,a) are dimensionless 
couplings. 
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In the physical spacetime dimension d = d + d (the continuation to complex dimensions will 
be discussed later) the classical theories with lagrangians are invariant under the weighted 

dilatation 

x^xe"^, x^xe"^/", ^ ^ ^ ^^i'^/'^-^) ^ (2.10) 

where O is a constant parameter. Each lagrangian term scales with the factor d, compensated by 
the scaling factor of the integration measure d'^x of the action. 

We call the models (|2.9p homogeneous. Homogeneity is preserved by renormalization, namely 
there exists a subtraction scheme in which no lagrangian terms of weighted degrees smaller than 
d{N) are turned on by renormalization. This fact is evident using the dimensional-regularization 
technique. Indeed, when (5^^ = d{N), the equality in (|2.8p holds, so io{G) = d{E) = 5^^\ 

The weighted scale invariance (|2.10|) is anomalous at the quantum level. The weighted trace 
anomaly and its relation with the renormalization group are studied in section 6. 

Nonhomogeneous theories are those that contain both weighted marginal and weighted rel- 
evant vertices. In these cases the weighted dilatation (|2.10p is explicitly broken by the super- 
renormalizable vertices, and dynamically broken by the anomaly. 

Let us analyze some explicit examples, starting from the homogeneous models. 

Homogeneous models. We begin with the (/?^-theories. Setting A'max = 4 in (|2.7|) we get 

10 

y<^<4. (2.11) 

One solution with d= 4 is the usual Lorentz-invariant (/9^-theory in four dimensions {d = d = 2, 
n = 1). A simple Lorentz- violating solution is the model with n = 2 described by the lagrangian 

^(2,4) = \{d^? + ^(AV^)' + (2-12) 

in six dimensions, with d = 2, d = 4. The (/7^-theories with n = 2 are used to describe the critical 



behavior at Lifshitz points [TTl [12] . 

It is clear that (12. lip admits infinitely many solutions for each value of d. For example, given 
a solution, such as (I2.12p . infinitely many others are obtained multiplying d and n by a common 
integer factor. For d= 4 we have the family of 2{n + l)-dimensional theories 

k2,2n) = + + ijXf^^'- (2.13) 

In general, for every Lorentz-invariant renormalizable theory there exists an infinite family of 
Lorentz- violating renormalizable theories. 
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Let us now concentrate on four dimensions. The spacetime manifold can be split as (d, d) = 
(0,4), (1,3), (2,2), (3,1), (4,0). There is no nontrivial solution with d = 0. Indeed, (12. Th implies 

4 



2 - n 

so n can only be 1, which gives back the Lorentz invariant V9^-theory. For d = 1 we get 



2 (n + 3) 
3 — n 



The only nontrivial solution is n = 2, which implies N^^ 



10 and 



(1,3) 2^ ^> ^ 2kV 6!Ar ^ ^' 10!A6 



-ip 



(2.14) 



For d = 2 we get N^a,x = 2 (n + 1): every integer n > 1 defines a nontrivial solution in this case. 
The simplest example is (d, d) = (2,2), n = 2. Listing all allowed vertices we get the theory 



^(2,2) = 2(5^) +^ 



4!Ai 



6!Ai 



(2.15) 



This model belongs to a family of d= 3, (2 + n)-dimensional (/3®-theories, whose lagrangian is 



C 



1 



1 



+ 



A. 



6!A 



2{n-l) 



(2.16) 



' if) 



6!A 



2(n-l) 



(2.17) 



when n is odd and 

when n is even. Observe that (I2.16P includes the Lorentz-invariant (/9^-theory in three spacetime 
dimensions, which is the case n = 1. 
For d = 3 we get 



2 (3n + 1) 



n + 1 

The solution with n = 2 has A'max = 4. However, this solution is trivial, since its unique vertex 
would have just one 9-derivative. Instead, for every n > 3, A^'max is equal to 5. For example, the 
theory with n = 3 is 

which is clearly unstable. Imposing the symmetry (/? — > —93 we have the modified 99^-theory 



1 



1 



A, 



which is stable for A4 > 0. Finally, for d = 4 we get again the Lorentz-invariant v9^-theory. 
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Nonhomogeneous models. Nonhomogeneous theories can be obtained from the homogeneous 
ones adding all super-renormalizable terms, which are those that satisfy the strict inequality 
< d{N). For example, keeping the symmetry — > —ip, the nonhomogeneous extension of 
(f2l2]l is just 

and the one of (I2.15P is 



Splitting the spacetime manifold into the product of more submanifolds. Instead of 
splitting the spacetime manifold into two submanifolds, we can split it into the product of more 
submanifolds, eventually one for each coordinate. This analysis covers the most general case. We 
still need to distinguish a d-dimensional submanifold containing time from the dj- dimensional 
space submanifolds , i = 1, . . . ^, so we write 

e 

4 = 1 

Denote the space derivatives in the ith space subsector with di and assume that they have weights 
1/nj. Then the kinetic term of the lagrangian reads 

I ^ 1 — 

where ^2(^1, A^) is the most general weighted homogeneous polynomial of degree 2 in the spatial 
derivatives, P2{X^^^^di, Al) = X'^P2{di,AL), invariant under rotations in the subspaces . The 
A^-dependence is arranged so that P2 has dimensionality 2. The previous analysis can be repeated 
straightforwardly. It is easy to verify that the weighted power-counting criterion works as before 
with 

Edge renormalizability. By edge renormalizable theories we mean theories where renormal- 
ization preserves the derivative structure of the lagrangian, but the powers of the fields are un- 
restricted. With scalars and fermions, such theories contain arbitrary functions of the fields and 
therefore infinitely many independent couplings. The notion of edge renormalizability is interest- 
ing in the perspective to study gravity. Indeed, Einstein gravity is an example of theory where all 
vertices have the same number of derivatives, but are nonpolynomial in the fluctuation around 



9 



flat space. Yet, diffeomorphism invariance ensures that the number of invariants with a given 
dimensionality in units of mass is finite. Therefore, in quantum gravity a polynomial derivative 
structure is sufficient to reduce the arbitrariness to a finite set of independent couplings. 

Edge renormalizable theories are those where u!{G) does not decrease when E increases, rather 
it is independent of E. By formula (12. 8p this means d= 2 {Nj^ax = oo), in which case uj{G) is 
always equal to 2. Since d= 2, d can be either or 1. The theories with d = contain higher 
time derivatives, so they are not unitary. Thus we must take d = 1. The homogeneous theory in 
four dimensions has lagrangian 

£ = £free + A, (2.18) 

where 

L 

and 

+V5{^)E\{di^f+ V(^{^)(Kipf + v-j{^){di'K^f + Vf,{ip){d^d,dkvf + v^9(v')aV, (2.19) 

where the V^'s are unspecified functions of ip with Vi{ip) = 0{ip)^ V4,{Lp),V-j{(p),V%{ip)^VQ{ip) = 
The lagrangian of the most general nonhomogeneous theory is (|2.18p with 

Aree = ^{dlff - -ip ( + ^ j 

and Ci equal to (|2.19|) plus 

^io(v') + Vii{p)Ep + yi2(v^)AV + vMiKpf + Vu{p)[{dipff, 

with Vii{p),Vi2{p) = 0{p^), Vvi{p) = 0{p). 
3 Inclusion of fermions 

In this section we classify the models of interacting fermions and scalars. We start from pure 
fermionic theories, with quadratic lagrangian 

1 „ 

£free = + TT^rrrV'^ ^, 

where n is the maximal number of 9-derivatives. The propagator 
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is, in momentum space, a weighted function of degree —1. The loop-integral measure is, as usual, 
a weighted measure of degree d. For the purposes of renormalization, the kinetic terms with fewer 
than n 5-derivatives can be treated as vertices. 

Label the vertices that have 2N ip-ip -legs by means of indices a and denote their weighted 
degree with 5^^\ Consider a diagram G with 2E external V'-V' -l^gs, constructed with -y^^ vertices 
of type {N,a). Once the subdivergences have been subtracted away, its overall divergence is a 
weighted polynomial of degree 



Lo{G) = d- 



E{d-l)+J2^P 

{N,a) 



d{l-N)-N 



in the external momenta. Renormalizability demands 

?(«) 



S'n' < d{l-N) + N = dF{N). 



(3.1) 



Polynomiality demands 

d>l, 

in which case the maximal number of external i/j-iIj -legs is 

N - ' 

-i'ma.v — 



d-1 



Pure fermionic homogeneous models have strictly renormalizable vertices, namely those with 
dpiN). Their lagrangian has the form 



A 



{N,a) 



C = ipBtp + ip+ / — ^ 



Here 



denotes a basis of lagrangian terms constructed with N fields ?/', N fields 
■0 and ndp{N) (^-derivatives, invariant under the reduced Lorentz symmetry. For simplicity, we 
can assume also invariance under parities in both portions of spacetime. 

Let us concentrate on four spacetime dimensions. The Lorentz split (1,3) gives A^max = 
1 -|- [n/3], which admits infinitely many nontrivial solutions, beginning from n = 3. For example, 
the n = 3 and n = 6 theories read 



Li (y L 



a ^ 



respectively. The Lorentz splits (2, 2) and (3, 1) do not admit nontrivial solutions, since A^max = 1 
in those cases. 
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Now we study the models containing coupled scalars and fermions. It is important to note that 
when different types of fields are involved, they must have the same n. We classify the vertices 
with labels {N^, N^, a), where 2N^ is the number of ^p-^p -legs, is the number of (/9-legs and a 
is an extra label that distinguishes vertices with different structures. Call ^ ^ the weighted 
degree of the a-th vertex. Consider a diagram G with 2E^ external V'-V' -legs, external f-legs 



and V 



(a) 



N ) vertices of type {N^,Nip,a). Once the subdivergences have been subtracted away, 
the overall divergent part of G a is a weighted polynomial of degree 



uj{G) 



d- 
+ 



E^{d-l)-^{d-2) 



E 



(a) 



in the external momenta. Renormalizability demands 



1 - N,, 



< 



1 - K 



+ N^ + N^ = d{N^,N^). 



Because 5^ 
inequality 



(«) 

{N^,N^) 



is nonnegative, the numbers of fermionic and bosonic legs are bound by the 



NJd- 



l) + ^{d-2)<d. 



Polynomiality demands d> 2. 

The homogeneous models have a lagrangian of the form 



A 



n-l 



1 ,^ 



2A 



2n-2 



(8V)- 



In four dimensions the splitting (1,3) has a unique nontrivial solution, which is the model (|2.14j) 
coupled to fermions. It has n = 2 and its lagrangian reads 

^(1,3) + ^MV- + \{M + ^(^v^)' + ^v^'OTV') + • (^TV-) 



'^4 AT i ^ 

+ .. . c, WW + 



4!Ai 



6!Ai 



Ai 



10!A| 
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The splitting (2, 2) admits infinitely many solutions. The simplest one is the theory with 
n = 2, symmetric under <-> —if, that couples (|2.15|) to fermions: 



'(2,2) 



The splitting (3, 1) admits, again, infinitely many solutions. 



■^4 ^2 I ^^6 6 

(f [oip) + 



4!Ai 



6!Ai 
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4 Renormalization 



In this section we study the structure of Feynman diagrams, their divergences and sub divergences, 
and the locahty of counter terms. For definiteness, we work with scalar fields, but the conclusions 
are general. 

One-loop. Consider the most general one-loop Feynman diagram G, with E external legs, / 
internal legs and f vertices of type {N, a) and weighted degree (5^^ . Collectively denote the 
external momenta by k. The divergent part of G can be calculated expanding the integral in 
powers of k. We obtain a linear combination of contributions of the form 

^ii'"Liji-j2.'^''i ■■■^'^uK--- K : (4.1) 



where 

To avoid infrared problems we insert a mass m in the denominators. For the purposes of renor- 
malization, it is not necessary to think of m as the real mass. It can be considered as a fictitious 
parameter, introduced to calculate the divergent part of the integral and set to zero afterwards. 
The real mass, as well as the other parameters Um of (12. 2p . can be treated perturbatively, so they 
are included in the set of "vertices". 

From the weighted power-counting analysis of section 2 we know that the numerator of (14. ip . 
namely 

Pij.1 ■ ■ ■ P^l2r Pji ' ' ' Pj2s ^I'l ' ' ' ^'^u ki^ ■ ■ ■ ki^ , 
is a weighted monomial Pq^n{p,k;p,k) of weight 1/n and degree 

n n ^-^ 

(N,a) 

At one loop the number of vertices equals the number of propagators. Using (12.3(1 and (5^^ < d{N) 



we get 



By symmetric integration, we can write 



r(^.") /v(2) j{I,n) ■j{I,n) 
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where S^i---i_i2r "^jr- j^s appropriately normalized completely symmetric tensors constructed 
with the Kronecker tensors of and M^, respectively. Performing the change of variables 

^2 X (n-l)/(2n) 



Pi=V^ri^] , (4.3) 



the integral Trlf'^ can be calculated using the standard formulas of the dimensional-regularization 
technique. We obtain 



n(47r)^/2r(D/2)r {D/2) T (/) 

The factor 1/n is due to the Jacobian determinant of the transformation (|4.3I) . The singularities 
occur for 

/<r+- + |. (4.4) 
n 1 

Combining this inequality with (|4.2I) we find that the divergent contributions satisfy 

u^'^<d^E[\-^=d{E). (4.5) 

The counterterms are a Pu+v/n,nif^j k): 

1- - - - £2 

- fevi • • • ki,--- h , where £ = d-D = ei-\ . 

£ n 

Thus (14. Sp ensures that the divergent terms can be subtracted away renormalizing the fields 
and couplings of the initial lagrangian. Observe that while the poles are proportional to 1/e, the 
residues of the poles can depend on ei and £2 separately. We know that taking a sufficient number 
of derivatives with respect to the masses, the external momenta and the parameters Om of (12. 2p . 
the integral becomes convergent. Therefore, the finite parts are regular in the limits £1,62 0, 
which can be safely taken in any preferred order. Objects such as ei/e and £2/£ multiply only local 
terms, so they parametrize different scheme choices and never enter the physical quantities. These 
observations generalize immediately to all orders. We define the minimal subtraction schemes as 
the schemes where 

£1 = a£, £2 = n{l — a)£, 
with a=constant, and only the pure poles in £ are subtracted away, with no finite contributions. 
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Overall divergences and subdivergences. Before considering Lorentz violating theories to 
all orders in the loop expansion it is convenient to briefly review the usual classification of diver- 
gences and the proof of locality of counterterms |17] in Lorentz symmetric theories. Consider the 
L-loop integral 

with Lorentz invariant propagators + m^), where k denotes the external momenta. The 

ultraviolet behavior of 2{k) is studied letting any (sub)set of the momenta p^^\ . . . ,p^^^ tend to 
infinity with the same velocity. Proper subsets of the momenta test the presence of subdivergences, 
while the whole set tests the presence of overall divergences, i) When any sub convergence fails, 
counterterms corresponding to the divergent subdiagrams have to be included to subtract the 
subdivergences. ii) Once all subdivergences are removed, the subtracted integral Zsub(^) can 
still be overall divergent. Taking an appropriate number M of derivatives with respect to the 
external momenta k the integral d^Isuh{k) becomes overall convergent. This proves the locality 
of counterterms. 

The overlapping divergences can be tested sending momenta to infinity with different velocities. 
For example, rescale pi, . . . ,pL as Xpi, . . . , Xpi, A^p;+i, . . . , X'^pi- This test, however, is already 
covered by the previous ones, since there is always a (sub)set Sfast of momenta tending to infinity 
with maximal velocity. In the example just given, Sfast = {Pi+i} ■ ■ ■ jPl)- The other momenta 
^siow grow slower, so they can be considered fixed in the first analysis and taken to infinity at a 
second stage. Weinberg's theorem [18] ensures that when Sfast tends to infinity the behavior of 
the relevant subintegral is governed by power counting and can generate logarithmic corrections 
depending on the momenta of Sgiow Then, when Sgiow tends to infinity the behavior of the integral 
over Ssiow is still governed by power counting, because the corrections due to the integrals over 
Sfast do not affect the powers of the momenta Sgiow Thus the power-counting analysis done in 
steps i) and ii) suffices. 

Now we generalize the analysis to Lorentz violating theories. We say that the components p 
and p of each momentum are rescaled with the same "weighted velocity" when 

p — > Xp, p A^/"p. 

Step i) is modified studying the convergence when any subset of momenta tend to infinity with the 
same weighted velocity. Whenever a sub convergence fails the counterterms associated with the 
divergent subdiagrams have to be included. Once the subdivergences are subtracted away, step 
ii) consists of taking an appropriate number of "weighted derivatives" (see below) with respect to 
the external momenta, to eliminate the overall divergences. It is easy to check that this procedure 
automatically takes care of the overlapping divergences. 
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Weighted Taylor expansion. Every Taylor expansion 



u\v\ 



ky^ ki^ ■ ■ ■ ki^ 



u=0 v=0 

can be rearranged into a "weighted Taylor expansion" 



f{Kk) = Y.j,f^'^{k,k\ 



where 



[i/n] 

= E „,„_„„),. 

is a weighted liomogeneous polynomial of degree ijn: 

f^'\Xk,X^/^k) = X'/y('\k,k). 

The £-th weighted derivatives with weight 1/n are the coefficients fui-uuM-ie-nu- 

The weighted Taylor expansion is useful to subtract the overall divergences. The overall- 
subtracted version of an integral whose weighted degree of divergence is lo reads 



(27r)^5 (27r)^^ 



nui ^ 

Q{p,p-,k,k) -J2-^Q^'Hp,p-Xk) 



where Q^^^ denotes the ^-th homogeneous polynomial of the weighted Taylor expansion of Q in 



Subtraction algorithm. Consider an L-loop diagram with V vertices and / propagators. The 
integrand, which we denote with Qg^ is a ratio of weighted polynomials and has degree equal to 
dq = X](jva) '^Tv'*'"?^^ ~ "^^^ integral X is a weighted function of degree dx = dg+DL. It has 
the form ^ _ 

d^^p f d^^p 



^Qg{p,P, k), 



(4.6) 



where p and p collectively denote the components of the momenta circulating in the loops, while 
k = {k, k) collectively denotes the external momenta. The overall degree of divergence of I is 
u{G) = dq+dL. 

The subtraction of divergences can be arranged according to the following table: 



QG{p,p;k,k) -T,jerQjiP^P'^k,k) 
~^t=o ^ TxQ^s {PtP'tk^k) S^=o ^ ir X^7er ^' ^) 



(4.7) 
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Here F denotes the set of divergent subdiagrams 7 of the diagram G. The rational function is 
obtained replacing the subintegrand with the appropriate, truncated, weighted Taylor expansion 
in the external momenta of 7. In the arrangement of (|4.7I) subdivergences are subtracted row-wise. 
Overall divergences are subtracted column-wise. 

A potential caveat comes from certain "extra subdivergences",those that occur when a subdi- 
agram 7' is convergent in Qg, but becomes divergent in one of the Q^'s. Then 7' does not belong 
to r, so its subdivergence is not subtracted row-wise. Nevertheless, it is easy to show that the 
extra subdivergences are automatically subtracted column- wise in (|4.7p . Details and an explicit 
example are given in appendix A. 

Thus, once the subdivergences have been subtracted away, the divergent part of every Feynman 
diagram is a weighted polynomial of degree (second row of (14.71) ) and can be removed 

renormalizing the lagrangian (|2.18|1 . 

5 Renormalization structure and renormalization group 

In this section we study the renormalization group. We illustrate it first in the d= 4 models 
(|2.13|) . For the reasons that we explain below, it is convenient to parametrize the bare lagrangian 
as 

^(2,2n)B = l^id^B? + + ^^^(„_i)(2_e./n) '^B (5-1) 

with 

ipB = Zl'''ip, Alb = ZaAl, Ab = A/u"Za, e = ei + -. (5.2) 

Observe that D= 4 — e. The weighted scale invariance (|2.10|) can be extended to a transformation 
that rescales also fi: 

x^xe-^, x^xe-^/", <^ ^ <^ e^(^-2)/^ // ^ /xe^. (5.3) 

The invariance under this transformation is not a symmetry. It just tells us that at the quantum 
level the weighted scale invariance (I2.10p is equivalent to a ^-rescaling. What is important in 
(|2.10l) and (jS.Sp is that A^, is unmodified. Because of (|5.3p . every renormalization constant in 
(|5.2p is just a function of A (otherwise it could also depend on evanescent powers of the ratio 
h/Al)- Thus, in the minimal subtraction scheme the A-beta function has the usual form 

The finiteness of (3x proves that all poles contained in Zx are inverse powers of e. 

In more detail, let us consider the contribution of a graph G with E external legs, / propa- 
gators and V vertices to the generating functional of one-particle irreducible diagrams. Such a 
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contribution has the schematic form 

7 ^V'(n-l)(2-e2/n) ^ 

where G denotes the value of the Green function. The dimensionahty of G in units of mass is 

[G]=d{v-^ + 1^ +E-4V, 

while its weighted degree is 

uj{G) = [G] - 6[G] =4-E + Auj{G), 

where 

Au;{G) = -e - I + , 5[G] = (2 - ^) (n - 1) - | + l) . 

Recalling that I is invariant under the weighted scale transformation (|5.3I1 , we find that G trans- 
forms as 

G ^ e^'^(^)G. (5.4) 

Once the subdivergences have been inductively subtracted away, the divergent part Gdiv is a 
weighted polynomial of degree 4 — in the external momenta. Matching the dimensionality and 
the weighted rescaling (15.41) we find 

where P4-E,n{d, 9; A^,) is a homogeneous weighted polynomial of degree A — E and dimensionality 
equal to its degree. The corresponding lagrangian counterterm reads 



where [P] means that the derivatives contained in P act on the scalar legs tp^ as appropriate. In 
particular, summing up all contributions for = 4, we get 



- [d^x , , y y^cL\\ 

'^L L=l 

where cl are divergent constants. Thus the renormalization constant of A is a power series in A, 

00 



L=l 
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with no spurious dependence on /u/A^. The same conclusion holds for the other renormalization 
constants. We have 

dAL . dhi^A 
= = -^h^- 

The Callan-Symanzik equation has the same form as usual. Calling 

Gk{xi, ■■■ • • • ,xk; = ((^(xi) • • • f{xk)) , 

we have 

d ^ d d \ 
fJ'-gj^ + ^^ax '^^^'^J ^k{xi,--- ,xk;xi,--- ,Xk;\,AL,iJ.) = 0. (5.5) 

The equation can be immediately integrated to give 

Gk{xi, ■ ■ ■ ,xk;xi, ■ ■ ■ ,xk; A, Al,^/^) = z^''{t)Gk{xi, ■■■ • • • ,xk; A(t), AL(t), /i), 

where t = In ^ and 

z{t) = exp (^*7^(A(t'))dt') , ^ = -Am), Aiit) = A^exp ^\^(A(t'))dt' 

with A(0) = A. Now the renormalization-group flow specifies how the correlation functions changes 
under a weighted overall rescaling. Indeed, the weighted scale invariance (|5.3|) - (|5.4|) tells us that 

Gk{xi,--- ,Xk;xi,--- ,Xfc;A,AL,e/i) =e^^'^"^^/^Gfc(exi,--- , ^Xfe; ^^/"xi, • • • , ^^/"xfc; A, A^, /i). 

A one-loop calculation for the models (I2.13P gives 

so these models are IR free. Only the beta function has a nonvanishing one-loop contribution. 
Indeed, using the dimensional-regularization technique tadpoles vanish in homogeneous models, 
so and rji start from two loops. 

Let us now consider the model (|2.15ll . The bare lagrangian reads 

^(2,2)B = \{dipBf + (Avpb)' + — ^§^</^|(5(/jb)' + -z^^A, 

where 

(/PB = ^yV, Alb = ^aAl, A4B = /x' (A4 + A4) , Agb = Ai^' (Ae + Ae) , e = ei + y. 

The theory is invariant under the scale transformation (|5.3I1 with n = 2. At one-loop we find 
Z^p = 1, = 1 and 

^ _ 5A4 ^ _ 5A4A6 5A4 



2(127r)2e' ° (87r)2e (487r)2£' 
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so the beta functions read 

n - r\ ^-^4 o _ , 5A4A6 5A| 

- + 2(12^' - + W? " (48^' 

The asymptotic solutions of the RG flow equations are 

2(127r)2 1 2 

A4~^^, A6~-A4, 

where t = In \x\^ and |x| is a typical weighted scale of the process. Since A4 and Ag must be 
nonnegative, the theory is IR free. 

6 Weighted trace anomaly 

The weighted scale invariance (I2.10p of the homogeneous models can be anomalous due to the 
radiative corrections. In this section we calculate the weighted trace anomaly, following [19]. For 
definiteness, we work with the model (|2.12l) . but the discussion generalizes immediately to the 
other models. 

Weighted dilatation. In the case of the model (|2.12|) . write the lagrangian as C{ip^ d^ip, Aip). 
The infinitesimal version of the transformation ()2.1Q|) reads 

6ip = Q ^ + X ■ d + ■ (f = QDip, 
with 0, <^ 1. The conserved Noether current = {J^, J^) is given by 

We continue the spacetime dimensions to complex values as explained in section 1. The continued 
transformation 5ip' and the continued current J'^ are obtained replacing Dip in 5ip and with 

fj^p = (^-l + x-d + ]^x-^Lp (6.1) 

(see (15. 3|] ). where D= 4 — e. At the bare level, the anomaly of (|6.ip is expressed by the divergence 
of J'^'. We find 

= -e^. (6.2) 
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Improved energy-momentum tensor and its weighted trace. The anomaly of the weighted 
dilatation is encoded also in the energy-momentum tensor, precisely in its "weighted trace". Let 
us start from the energy-momentum tensor given by the Noether method. For the model (|2.12l) 
we have 

dC dC ^ 

d{d^ip) d{Aip) 

This tensor is not symmetric, but conserved: it is easy to check that d^T^^u = 0, using the field 
equations. Next, define the improved energy-momentum tensor 

^ ^ 1 . a r - 2 ^ 3D-2DD + 3D-5 _ , - . 

Tf,u = o^tpd^ip - -^dytp df,Aip - df,uC - — ^ — —^t^uV + 7^ — ttt^ -^tiu [V^V) 

4(D - I) [-U - ijA^ 



3 -2D _ ,^ .2 3 -2D 



2"i^u{daV) H T^tiai^T^au^p) ■ (6.4) 



2 

where 7?^,^ = 9^9^ — dfj_,^d^ and W^^ = d^d^ — 6fj,,yd . The first three terms of (|6.4p correspond 
to the Noether tensor (|6.3|1 . while the rest collects the improvement terms, identically conserved. 
Define the weighted trace 

Using the field equations, it is easy to show that Tf^u is conserved and that its weighted trace 
O vanishes in the physical spacetime dimension d = d + d. Moreover, T^y is conserved also in 
the continued spacetime dimension. The coefficients of the improvement terms are chosen so 
that in the free-field limit vanishes also in the continued dimension D = D + D. Finally, it 
is straightforward to check that the weighted trace coincides with the divergence (|6.2p of the 
current J'^. 



Anomaly. We need to write in terms of renormalized operators. When we differentiate a 
renormalized correlation function with respect to A or A/, we obtain a renormalized correlation 
function containing additional insertions of —dS/dX or -dS/dAi, respectively. Thus, —dS/dX 
and -dS/dh-L are renormalized operators. Following a standard procedure [H] we can find which 
operators O they are the renormalized versions of. In the minimal subtraction scheme, it is 
sufficient to express the renormalized operators as bare operators Ob plus poles. Schematically, 

finite = Ob + poles =^ finite = [O] . 

where [O] denotes the renormalized version of the operator O. We find 

dS , . I f ,^ , , dS Xb f A fi' f, 4, 
dX = ^"^'^ = % V^-f^-] - ^^^^MZ - "4!A|7 y ^^J = 4!Ai y 
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where [E,^] = j ip{6S/5ip) is the operator that counts the number of (/^-insertions. Thus, 

y ® = -y ^4!A|^ = 4!Ai y[H-^y[(A^)]-7.[i?.]. 

The result agrees with the Callan-Symanzik equation (|5.5p . which can be expressed as 

e (^(xi) • • • fixk)^ = ivixi) ■ ■ ■ fixk)) . 

Indeed, 

dS dS ^ dS 



7 Nonrelativistic theories 

Nonrelativistic theories can be studied along the same lines. The action contains only a single 
time derivative d, 

SO the theory is more divergent. The dimensional-regularization is not easy to use, since there 
is no simple way to continue the single-derivative term '(fdf to complex dimensions. Thus we 
assume an ordinary cut-off regularization. 

The propagator is defined by the term Tpdf plus the lagrangian quadratic term with the highest 
number of 9-derivatives, say 2n, 

Ciree=^(^d+^^^ if. 

For the purposes of renormalization, the other quadratic terms, if present, can be treated per- 
turbatively, as explained in section 2. Thus the nonrelativistic propagator is the inverse of a ho- 
mogeneous weighted polynomial of degree 1 and weight 1/n. The integral measure has weighted 
degree d= 1 + (d— l)/(2n). A Feynman diagram G with E total external legs, I propagators and 
v^^^ N-leg vertices of weighted degrees is a weighted function of degree 



iN,a) 



Formulas (12.31) still hold. We have 



ujiG) = d-^id-l)+ 

(N,a) 



(a) f N \ N' 



(a) 
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Renormalizable theories are those that contain the vertices with 

N fN 

Strictly renormahzable theories are those that have 



4"^<^-(^-l)'^- (7.1) 



jI?) = 1, - ( - 1 u. 



N fN 



^ 2 V 2 
Polynomiahty requires now 

which ensures also that decreases when the number of external legs increases. The maximal 
number of legs is 



-^max 

It is straightforward to check that E = N implies 



2d 
d-l 



(7.2) 



u{G)<d-^{d-l), 

so by (|7.ip the type of vertex that subtracts the divergence of G is already present in the la- 
grangian, which proves renormalizability. No terms with more than one time derivative are turned 
on by renormalization. 

Let us now see some examples of homogeneous models, beginning from the (/?^-theories. Setting 
A^niax = 4 in dlJ) we get 

5<<i<2. 

For d= 2 we have d = 2n + 1 and the family of odd-dimensional theories 

Setting A^max = 6 we have 7/5 <d< 3/2. For d= 3/2 we have d = n + 1. If n is odd we have the 
family 

In particular, we see that there exist four-dimensional (n = 3) nonrelativistic renormalizable 
99^-theories. If n is even we must include additional vertices, 
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8 Conclusions 



In this paper we have classified the unitary Lorentz violating renormalizable quantum field theories 
that can be obtained improving the UV behavior of propagators with the help of higher space 
derivatives. The removal of divergences is governed by a weighted power-counting criterion. If 
the lagrangian has an appropriate form, time derivatives are "protected", in the sense that no 
higher time derivatives are turned on by renormalization. The so-defined theories are unitarity, 
but have modified dispersion relations. We have studied their main properties, including the 
renormalization group flow and the weighted trace anomaly. 

Natural extensions of this work are those that aim to include gauge fields and gravity. Pos- 
sible applications range from high-energy physics, effective field theory, nuclear physics and the 
theory of critical phenomena. In the high-energy physics domain, it would be interesting to 
explore the work-hypothesis that Lorentz invariance is violated at very high energies, to define 
the ultraviolet limit of quantum gravity, or study new types of Lorentz invariant extensions of 
the Standard Model. It would also be interesting to embed the weighted scale invariance into a 
"weighted conformal group", generalizing the Galilean conformal group that characterizes a class 
of nonrelativistic theories |20| . 



Appendix A: extra subdivergences 

In this appendix we give more details on the extra subdivergences mentioned in section 4. By 
construction, every row of table (|4.7p is free of "ordinary" subdivergences, namely those originated 
by the subdiagrams 7. Every column is free of overall divergences. Extra subdivergences are 
those that occur when a subdiagram 7' is convergent in Qg, but becomes divergent in one of the 
(5^'s, after replacing 7 with its counterterms. Here we prove that the extra subdivergences are 
automatically subtracted column-wise. 

It is useful to have an explicit example in mind, such as the two-loop diagram depicted in fig. 
1, in the four dimensional (/^'^-theory. The diagram is the p-q integral of 

^ 1 1 

{p^ + m?) [{p — k)^ + m?] [q^ + m?) [{q+ p + k')"^ + m?] 
The p-integral is convergent, the g-integral is not. The sub divergence is subtracted by 

{p^ + 'm?) [{p — k)"^ + m?] {q^ + m?)"^ ^ ^ 

In this expression, however, the p-integral is divergent. This divergence is what we call an extra 
subdivergence. The table reads 
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Figure 1: Simple example of diagram that generates an extra sub divergence 



In the general case, assume that the subdiagram 7 contains I loops and that contains 
some extra subdivergences. The extra subdivergences can be overall or not. We call them overall 
if they arise letting all of the remaining L — I loop momenta tend to infinity. They are not 
overall if they arise letting only a subset of the remaining L — I loop momenta tend to infinity. 
Proceeding inductively, we can assume that the non-overall extra subdivergences have already 
been subtracted away. Thus, we need to consider only the overall extra subdivergences. It is not 
difficult to see that they are subtracted column- wise in (|4.7p . Indeed, as in (jA.lj) . the integrands 
that generate extra overall subdivergences factorize (or split into a sum of terms each of which 
factorizes): one factor is responsible for the extra subdivergence (see the first factor of —Q^ in 
(jA.ll) ). while the other factor is the 7-counterterm (see the second factor of —Q^ in (jA.ip ). The 
second factor is the same throughout the column. Thus, the column subtracts away the overall 
divergence of the first factor, which is precisely the extra subdivergence. Recapitulating, the rows 
are free of ordinary subdivergences and the columns are free of extra subdivergences and overall 
divergences. Thus the table (14.71) is convergent. In the example (IA.2j) . it is clear that the column 
of —Qj is p-convergent. 

Appendix B: Euclidean propagators 

Let us examine some propagators 



in coordinate space. The Euclidean (2,2)-propagator in four dimensions with n = 2 reads 



where I denotes the modified Bessel function of the first kind, while SL denotes the modified 
Struve function. For \x\ » A^^x^ and \x\ <^ Kijjp' we have 



G(2,2){x,X,Kl) 



Al 



^ [/o(AixV4!5?|) - 5Lo(AixV4|x|)] 



16|x 




Al 



and G(2,2) 



1 



16|a; 



27rx2' 



respectively. 
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Instead, the Euclidean (l,3)-propagator with n = 2 reads 



G(i,3)(2;,x, Al) = ^-p|Erf 



■2 



87r|x| \ y 4|ir[ / 
In the two hmits considered above we have the behaviors 

^3/2 ^ 

respectively. 
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